AD-A183  492 
UNCLASSIFIED 


THE  APPLICATION  OF  STATISTICAL  TURBULENCE  THEORY  TO  1/1 

CONVECTIVE  INSTABILITIES^)  SRI  INTERNATIONAL  HEHLO 
PARK  CA  C  L  RINO  81  NOV  86  DNA-TR-86-482 
DNA881-86-C-8882  F/G  4/1  NL 


AD- A 183  492 


OTIC  FILE  CORY 


DNA-TR-86-402 


THE  APPLICATION  OF  STATISTICAL  TURBULENCE  THEORY 
TO  CONVECTIVE  INSTABILITIES 


C.  L  Rino 

SRI  International 

333  Ravenswood  Avenue 

Menlo  Park,  CA  94025-3434 


1  November  1986 


Technical  Report 


CONTRACT  No.  DNA  001-86-C-0002 


Approved  for  public  release; 
distribution  is  unlimited. 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMC  CODE  B3220867663  RW  RA  00001  25904D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305-1000 


DISTRIBUTION  LIST  UPDATE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


1».  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2*.  SECURITY  CLASSIFICATION  AUTHORITY 

N/A  since  Unclassified 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

N/A  since  Unclassified 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERISI 

SRI  Project  1700 


6a  NAME  OF  PERFORMING  ORGANIZATION 

SRI  International 


6c  ADDRESS  /Cry.  Sum.  end  ZIP  Code I 

333  Ravenswood  Avenue 

Menlo  Park,  California  94025-3434 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


6b.  OFFICE  SYMBOL 
(if  applicable) 


3  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBERISI 

DNA-TR-86-402 


7.  NAME  OF  MONITORING  ORGANIZATION 

Director 

Defense  Nuclear  Agency 


7b  AODRESS  I City .  Stele,  end  ZIP  Code I 

Washington,  DC  20305-1000 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


8c  AODRESS  (City.  Stele,  end  ZIP  Code I 


8b  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(if  applicable) 


RAAE/Wittwer 


DNA  001-86-00002 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


project 

TASK 

NO 

NO 

RW 

RA 

62715H 


11  TITLE  (Include  Security  Classification) 

THE  APPLICATION  OF  STATISTICAL  TURBULENCE  THEORY  TO  CONVECTIVE  INSTABILITIES 


WORK  UNIT 
ACCESSION  NO 

DH200891 


12  PERSONAL  AUTHORIS) 

Rino,  Charles  L. 


13a  TYPE  OF  REPORT 

Technical 


13b  TIME  COVERED 

FROM  850313  TO  860930 


14  DATE  OF  REPORT  lYeer.  Month.  Deyl  15  PAGE  COUNT 

861101  32 


16  SUPPLEMENTARY  NOTATION 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMC  Code 
B3220867663  RW  RA  00001  25904D. 


17, _ COSATI  CODES _  18  SUBJECT  TERMS  (Continue  on  reverse  tf  necessary  and  identify  by  block  number ) 

field  group  sub-group  Turbulent  Plasmas 

T7  02  1  Convective  Plasma  Instabilities 

20 _ 09 _ Irregularity  Power  Spectral  Density _ 

J9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

:^kStatistical  turbulence  theories  have  been  used  effectively  to  characterize 
spectral  density  functions  in  fluid  turbulence,  and  in  turbulent  plasmas,  magnetic 
fields,  and  passive  scalars  in  turbulent  fluids.  Irf  t^his  report  review  the 
applications  to  convective  instabilities  of  the  turbulence  theories  ideve loped  original¬ 
ly  and  independently  by  Kraichnan  and  Kadmotsev,  bu^Pgeneral i zed  for  ionospheric 
appl icat ions , by  Sudani  These  instabilities  are  believed  to  be  the  dominant  structuring 
mechanism  in  equatorial  spread  F,  artificial  barium  clouds,  and  in  high-altitude 
nuclear  plumes  at  late  times.  We  hav^1  emphasi zed, the  ramifications  of  the  cross-field 
anisotropiy  unique  to  convective  instabilities,  which  have  not  been  treated  previously. 
We  show  that  the  limiting  form  of  the  one-dimensional  spectral  densitv  functions 
measured  perpendicular  to  the  F  x  B  drift  direction  is  proportional  to  k  where  k  is 
the  wave  number.  ,  • 


20  oistribution/availability  of  abstract 
Ounclassifieo/unlimiteoQsame  as  rpt 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Sandra  E.  Youn 


DO  FORM  1473,  84  mar  83  apr  < 


I  OTIC  users 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  tlnchnle  Ares  Codel\  22c  OFFICE  SYMBOL 


(202)  325-7042 


DNA/CST1 


S3  APR  edition  m*y  b«  ui #tJ  until 

All  Oth«r  •ditiont  »r«  ob«ol«t« 


SECURITY  CLASSIFICATION  of  this  page 


CONVERSION  TABLE 


Conversion  factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement 


i-  FY 

TO  GET 

MULI  Is  LT  P 

RY  ^ _ 

nivinr 

TO  GET  ^  . 

aagatrora 

1.  000  000  X  S  -10 

me  tan  (m) 

atmosphere  (normal) 

1  013  25  X  C  42 

kilo  pascal  (kPa) 

bar 

1  000  000  X  E  +2 

Idlo  pascal  (kPa) 

bam 

1  000  000  X  E  -28 

mater*  (m2) 

British  thermal  unit  (thermochemical) 

1.  OSS  350  X  E  43 

Joule  (J) 

calorie  (thermochemical) 

4. 104  000 

Joule  (J) 

cal  (thermochemical) /cm 

4. 104  000  X  E  -2 

mega  Joule/m2  (MJ/m2) 

curie 

3.  TOO  000  X  E  4l 

•gig  a  beoquerel  (GBq) 

degree  (angle) 

1  745  329  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

*  (ft  4  4  50  S7)A  8 

degree  ketvtn  (K) 

electron  volt 

1  802  10  X  E  -IS 

Joule  (J) 

ct% 

1.000  000  X  E  -T 

joule  (J) 

erg/aecood 

1.000  000  X  E  -1 

watt  (W) 

foot 

3.  048  000  X  E  -1 

meter  (m) 

foot -pound -force 

1.355  SIS 

Joule  (J) 

gallon  (l  S  liquid) 

3  785  412  X  E  -3 

meter2  (ro2) 

inch 

2  540  000  X  E  -2 

meter  (m) 

|erk 

1  000  000  X  E  4  9 

joule  (J) 

joule/ki lagram  <J/kgl  (radiation  doae 

absorbed) 

1 . 000  000 

Gray  (Gy) 

kilotons 

4  183 

lerajoules 

tup  (1000  Ihf) 

4  44  8  222  XE  43 

newton  (N) 

kip/inch2  (kai) 

6  804  757  X  E  43 

kilo  paacal  IkPa) 

ktap 

1  000  000  X  E  4  2 

newton  -second  /m 
(N-e/m2) 

micron 

1  000  000  X  E  -6 

meter  (m) 

mil 

2  540  000  X  E  -5 

meter  (m) 

mile  (international) 

1.  609  344  X  E  43 

meter  (m) 

ounce 

2  634  952  X  E  -2 

kilogram  (kg) 

pound -force  (lbs  avoirdupois) 

4  448  222 

newton  <N) 

pound -force  inch 

1  129  848  X  E  -1 

newton-meter  (N«m) 

pound  -force /uich 

1  751  268  X  E  42 

newton /meter  (N/m) 

pound  -force /toot2 

4  788  026  X  E  -2 

kilo  paacal  (kPa) 

pound  -force  /inch2  (psi) 

6  894  757 

kilo  paacal  (kPa) 

pound -mass  (Ibm  avoirdupois) 

4  535  924  X  E  -1 

kilogram  (kg) 

pound -mass -foot2  (moment  of  inertia) 

kilogram -meter 

4  214  011  XE  -2 

(l«-in2) 

pound -mass /Toot' 

kilogram /meter 

1  601  846  X  E  4  1 

(%/m3) 

rad  (radiation  doae  absorbed) 

1  000  000  X  E  -2 

••Cray  (Gy) 

roentgen 

coulomb^dl<gnm 

i 

2  579  760  X  E  -4 

(C/kg) 

shake 

1  000  000  X  E  -8  j 

second  (a) 

.  1 

slug 

1  459  390  X  E  «l 

kilogram  <%) 

torr  (mm  Hg,  0*  O  J 

1  .133  22  X  E  -1  j 

kilo  paacal  IkPt) 

•The  becquere!  (Bq)  ia  the  SI  unil  of  radioactivity,  1  Bq  1  event/s 
••The  Cray  (Gy)  is  the  SI  unit  at  absorbed  radiation 


SECTIOH  1 


INTRODUCTION 

The  simplest  mathematical  model  that  contains  the  essential  physics 
of  a  convectively  unstable  plasma  is  given  by  the  coupled,  nonlinear 
differential  equations 

JJ  +  (bxV*/B)  .  VN  =  KVN  ,  (1) 

and 


V(NV$)  =  D  .  VN  (2) 

where  N(p,t)  is  the  electron  density  in  a  plane  normal  to  the  magnetic 
field.  The  vector  D/B  is  the  the  cross-field  drift  velocity  caused  by 
the  combined  forces  of  gravity,  neutral  winds,  and  electric  fields. 

Thus,  D  can  be  thought  of  as  an  effective  electric  field.  The  poten¬ 
tial,  +(p,t),  which  is  established  by  Eq.  (2)  from  the  instantaneous 
electron  density  configuration,  modifies  the  effective  electric  field  in 
such  a  way  that,  as  time  evolves,  small-scale  structures  are  drawn  out 
of  the  trailing  edge  of  high-density  regions  and  are  thereby  enhanced. 
The  structure  is  ultimately  removed  by  cross-field  diffusion  at  a  rate 
determined  by  the  diffusion  coefficient  K. 

The  physics  of  this  two-dimensional  model  and  its  potential  limita¬ 
tions  are  discussed  in  detail  in  Perltlns  et  al.  [1973],  Ossahow  [1981], 
and  Zalesafc  et  al.  [1985].  Here,  we  shall  discuss  the  applications  of 
statistical  turbulence  theory  to  determine  the  quasi-steady-state 
spectral  density  function  (SDF)  that  characterizes  the  late-time 
evolution  of  plasma  structure  implied  by  Eqs.  (1)  and  (2).  As  yet,  no 
viable  theory  has  emerged  that  can  explain  the  measured  fcilometer-scale 


spectral  characteristics  of  naturally  occurring  equatorial  F-region 
ionospheric  irregularities  that  severely  affect  satellite  radio  signals 
r Aarons.  1982].  Yet,  it  is  generally  believed  that  the  development  and 
evolution  of  these  irregularities  can  be  modeled  by  using  straight¬ 
forward  modifications  of  Eqs.  (1)  and  (2)  fZalesak  et  al. .  1982].  Our 
motivation  in  undertaking  this  study  was  to  determine  if  formal  turbu¬ 
lence  theory  can  be  modified  to  explain  the  results  of  recent  data 
analyses . 

Ott  and  Parley  [1974]  first  pointed  out  that  Eq.  (1)  is 
mathematically  identical  to  the  two-dimensional  Navier-Stokes  (NS) 
equation  when  it  is  written  in  terms  of  vorticity, 

w  =  V  x  v  ,  (3) 

where  v  is  the  incompressible  fluid  velocity  (V.v  =0).  If  v  is  con¬ 
fined  to  a  plane,  it  is  easily  shown  that  the  momentum  balance  of  iner¬ 
tial,  pressure,  and  viscous  forces  can  be  written  as 

f*  +  v.Vw  =  vV2w  ,  (4) 

at 

where  v  is  the  kinematic  viscosity.  Thus  the  momentum  balance  for  two- 
dimensional  fluids  when  expressed  in  terms  of  scalar  vorticity  is  iden¬ 
tical  to  the  continuity  equation  for  electron  density;  however,  to 
complete  the  formal  parallel  we  must  identify  the  electron  drift 
velocity 

v  =  a  x  V4/B  (5) 

e  z 

with  the  fluid  velocity  v. 

Because  v  is  divergence  free,  it  too  can  be  derived  from  a  poten¬ 
tial  function.  If  we  let  v  =  a  x  74,  it  follows  from  Eq  (3)  that 

z 

=  -W  (  6 ) 


Thus,  for  inertially  driven,  two-dimensional  turbulence,  the  potential 
function  that  generates  the  velocity  structure  is  related  to  the  scalar 
vorticity  by  Poisson's  equation.  Because  of  the  profound  differences 
between  Eqs.  (2)  and  (6),  theories  that  address  the  turbulent  structure 
of  w  cannot  be  applied  directly  to  the  convective  plasma  system  de¬ 
scribed  by  Eqs.  (1)  and  (2).  Indeed,  even  if  one  were  to  accept 
ordinary  turbulence  as  the  source  of  the  structure  in  the  velocity 
field,  the  electron  density  must  act  as  a  passive  scalar  to  map  the 
inertially  driven  velocity  structure  rKellev  and  Ott.  1978].  Ideally, 
one  would  proceed  from  a  general  plasma  theory  that  accommodates  both 
collisional  and  inertial  forces  [Klntner  and  Sevier.  1985];  however,  the 
collisional  limit  is  the  most  difficult  to  treat,  and  it  appears  to  have 
some  unique  properties. 

Consider  that  in  Kolmogoroff 's  theory  of  fully  developed  turbu¬ 
lence,  locally  unstable  eddies  spontaneously  break,  up  into  smaller 
eddies  until  they  are  removed  by  viscous  dissipation.  In  sharp  con¬ 
trast,  convectively  unstable  density  irregularities  evidently  do  not 
spontaneously  break  up  or  bifurcate.  Rather,  early  in  the  nonlinear 
phase  of  the  structure  development,  fingerlike  strlations  develop  with  a 
comparatively  small-size  distribution.  Evidently  the  strlations  con¬ 
tinue  to  elongate  without  further  significant  change  in  their  individual 
shapes  r McDonald  et  al..  1981;  Zalesak.  et  al..  1985].  This  process 
has  been  referred  to  as  "freezing."  Clearly,  the  structure  dynamics  are 
vastly  different  in  the  collisional  and  inertial  models,  and  we  should 
expect  these  differences  to  manifest  themselves  in  a  viable  theory. 

Formal  plasma  turbulence  theories  have  been  developed  by  Kraichnan 
and  Kadmotsev.  Both  theories  use  the  Direct  Interaction  Approximation, 
which  was  developed  by  Kraichnan  [1958,  1959];  his  turbulence  theory  is 
usually  referred  to  as  the  DIA.  Kadmotsev* s  method,  which  was  developed 
more  Intuitively,  is  usually  referred  to  as  the  the  Weak  Coupling 
Approximation  (WCA)  [Kadmotsev.  1965].  Sudan  and  Pflrsch  [1984]  have 
shown  that  the  two  methods  ultimately  give  rise  to  formally  identical 
mathematical  models;  however,  the  WCA  form  suggests  a  somewhat  simpler 


procedure  for  deriving  spectral  characteristics.  Following  Sudan  and 
Pfirsch's  approach,  we  shall  start  with  the  DIA  because  it  proceeds 
directly  from  Eq.  (1).  We  shall  then  recast  the  equations  in  the  WCA 
form  and  follow  Sudan's  prescription  for  deriving  the  form  of  one- 
dimensional  spectral-density  function.  The  main  difference  is  that  we 
have  attempted  to  relax  the  assumption  of  isotropy,  which  is  critical  in 
all  the  applications  of  the  method  to  date  f Sudan  and  Keskinen.  1979; 
Sudan.  1983:  Sudan  and  Kesfclnen.  1984].  This  is  achieved  by  approxi¬ 
mating  the  anisotropy,  which  manifests  itself  in  the  "frozen"  strlation 
distribution,  by  a  one-dimensional  transfer  of  scalar  variance  in  the 
Fourier  domain.  This  provides  an  alternate  means  of  obtaining  a  one- 
dimensional  conservation  equation  in  place  of  the  angle-averaged  form 
that  is  normally  used. 

The  first  step  is  to  transform  Eq.  (1)  into  the  spatial  Fourier 
domain.  We  assume  that  M(p,t)  and  +(p,t)  admit  spectral  decompositions 


of  the  form 


N(p,t)j  |  N(k.t)  ^ 

j -  E<  ;exp  U*-p1 

♦(p.t)  \  /  iu,t)  \ 


where 


k  =  (nAK  ,  mAK  ) 
x  y 


for  all  Integer  values  of  n  and  m.  Thus,  N  and  4  are  periodic  on  the 

lattice  bounded  by  L  =  2i/AK  and  L  =  2i/AK  .  As  discussed  by 

x  x  y  y 

Kralchnan  [1975],  this  imposes  no  limitations  of  practical  significance. 
Also,  because  N  and  +  are  real,  N  and  4  have  the  hermitian  property 
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Substituting  Eq.  (7)  into  Eq.  (1)  and  performing  some  standard 
Fourier  manipulations  gives  the  spectral  domain  form  of  Eq.  (1),  namely. 

+  k2KN(k)  =  23  (b  x  k'  x  k"/B)  ♦(k,)H(kt ')  .  (10) 

O  t 

The  symbol  A  means  that  only  those  modes  that  satisfy  the  triangle 
equality 

k  -  k*  +  k* •  (11) 

contribute  to  the  summation.  We  could  have  included  the  delta  function 
5(k,  -  k*  -  k* ' )  in  the  summation,  but  the  alternate  notation  is  more 
compact.  In  the  interest  of  notational  efficiency,  we  have  also  sup¬ 
pressed  the  Independent  time  variable  t. 

A  * 

next,  we  multiply  Eq.  (10)  by  H  (k)  to  obtain 
•  2 

9ly(k)J_  +  2Kk2[M(k)|  =  E(k)  ,  (12) 

ot 

where 

A  . 

E(k)  =  2  23  <b.k'  X  k”/B)  (fr(k’)N(k”)  .  (13) 

Equations  (12)  and  (13)  characterize  the  temporal  evolution  of  the 
spectral  Intensity  for  a  single  realization  of  the  process.  The 
nonlinear  mixing  term  E(k)  has  the  important  property 

£E(k)  =  0  ,  (14) 

which  does  not  depend  on  how  $(k)  is  related  to  N(k).  Thus,  In  Eqs.  (1) 
and  (4)  the  nonlinear  interaction  of  the  Fourier  modes  does  not  change 
the  scalar  variance 


(15) 


Q  =  £|N(lt)‘ 


or  enstrophy  (the  name  given  to  the  variance  of  scalar  vorticity)  in  the 
NS  system.  This  conservation  property  is  the  basis  for  developing  a 
statistical  turbulence  theory. 

A  formal  statement  of  the  conservation  of  scalar  variance  can 
obtained  from  Eqs.  (12)  and  (14)  as 


+  2K  2>2|N(k)2|  =  0 


Thus,  the  structure  evolution  is  such  that  the  only  change  in  scalar 
variance  is  through  diffusion.  Turbulence  theory  attempts  to  charac¬ 
terize  the  average  spectral  characteristics  in  the  quasi  steady  state 
where  the  tendency  for  structure  to  form  is  essentially  balanced  by  the 
structure  removal  process.  In  fluid  turbulence,  energy  must  be  supplied 
continuously  to  sustain  this  steady  state.  In  collisional  plasma  struc¬ 
turing,  a  reservoir  of  structure  at  low  spatial  frequencies  can  sustain 
structure  growth.  Thus,  In  our  application  of  the  DIA,  we  do  not 
Include  a  forcing  function. 

We  envision  many  realizations  of  the  process  with  small,  randomly 
Irregular  changes  in  the  Initial  conditions,  but  with  similar  structure 
reservoirs  at  small  spatial  wave  numbers  to  sustain  the  growth  of  small- 
scale  structure  at  late  times.  We  expect  the  ensemble  average  of  the 
spectral  intensities  from  each  realization,  formally 


♦  (It)  =  <  |N(k)  |  >  ,  (17) 

to  yield  an  identifiable  SDF,  although  it  is  well  k.nown  that  Eqs.  (1) 
and  (2)  will  not  support  a  true  steady  state.  If  we  attempt  to  evaluate 
Eq.  (17)  by  using  Eqs.  (12)  and  (13),  we  encounter  the  third-order 

*  *  a  • 

moment  <<fr(k' )N(k' * )N  (fc)>.  Iterating  the  procedure  will  always  intro¬ 
duce  a  moment  of  higher  order.  Thus,  irrespective  of  the  mathematical 
subtleties  associated  with  statistical  stationarlty ,  the  formal  moment 


ww 
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equations  are  indeterminant.  Direct  interaction  is  a  closure  approxi 
mation  in  which  the  third-order  moments  are  evaluated  in  terms  of 
fourth-order  moments  containing  the  temporal  response  function  of  Eq. 
(1)  to  small  changes  in  $.  The  fourth-order  moments  are  evaluated  by 
using  the  well-known  relation  for  Gaussian  processed,  although  the 
Gaussian  assumption  itself  is  not  critical  to  the  development. 


SECTION  2 


♦  00  =  W(lO  N(lO  . 


( 


where  W(lt)  *  W  (—It.)  to  preserve  the  hermitian  property  of  ♦(*).  This  is 
tantamount  to  linearizing  the  potential  equation,  Eq.  (2),  but  within 
the  context  of  the  DIA,  this  is  a  local  approximation  that  preserves  the 
essential  nonlinearity  in  the  problem. 

* 

The  lowest  order  term  in  the  expansion  of  -N  is  assumed  a  Gaussian 
variable,  and  the  corresponding  term  in  the  expansion  of  G  is  assumed 
exact.  After  considerable  algebraic  manipulations,  we  obtain  the 
following  equations  for  the  time-dependent  SDF: 


,  as  (k;t,t') 

A  — fi -  +  —  ^ 

2  at 


- - +  Kk*sN(k,t,t')  =  £(b.k'  x  k'7B)2f(k',k") 

X  [j  G(k,t',t")  SH(k';t,t")  SH(lt";t.t")  dt" 


-  f 


G(*\t,t")  SN(k;t,t")  SN(k.";t',t")  dt"]  (22) 


where 


[ft  +  =  -  £(b.k'  X  k"/B)'f(k',k") 

|  G(k";t,t")  SN(k';t,t")  G(k;t'.t")  dt"  +  S ( t  -  t')  .  (24) 


f(k',k")  =  ^lW(k')  -  W(k")|2 


2 

**  w  *  •  /n  \  **  a 


Sjj(k.;t,t' )  =  <N(k;t)N  (k;t')> 


As  the  DIA  equations  are  written  here,  they  accommodate  both  the  NS  and 
convective  plasma  equations.  For  NS,  it  follows  from  Eq.  (6)  that 


-  *  V  V  ’J-  'J'V  •>  V  V  . 


W(k)  *  -k  ‘  .  (26) 

which  is  exact.  For  a  convectively  unstable  plasma,  it  follows  from  Eq. 
(2),  essentially  by  dropping  the  nonlinear  term,  that 


W(k)  =  -i  *5° 
It 


Because  Eq.  (26)  is  Isotropic,  the  functions  in  Eqs.  (22)  and  (24) 
depend  only  on  the  magnitude  of  It  and  the  equations  can  be  reduced  to  a 
one-dimensional  form  by  integrating  over  all  directions.  When  Eq.  (27) 
applies,  however,  the  basic  equations  admit  no  obvious  simplification. 

Thus,  we  transform  the  complete  DIA  equations  to  their  more  compact 
WCA  form.  We  first  assume  that  the  time  dependences  are  functions  only 
of  the  corresponding  time  differences,  whereby  Eqs.  (22),  (24),  and  (25) 
can  be  transformed  to  the  temporal  frequency  domain.  After  some 
straightforward  manipulations,  we  obtain  the  equivalent  WCA  equations 


2  °  «  p 

(w  -  iKk  )  SN(k;w)  =  -  £(b.k'  X  k"/B)  f(k';k") 


f" 

1  S  (k' ;w' )  S  (k' ' ;w  -  w' 
.  _ 25 _ n _  w _ 

[w  -  iKk2  +  r(k;w) ] 


(  S  (k";w  -  w')dw'  ) 

3  (k;w)  - “ - - -  < 

l  [w'  -  iKk'  +  r(k';w')]  ) 


“  2 

r(k; w)  =  £(b.k'  *  k"/B)  f(k';k") 


f  Sw(k";w  -  w')dw' 

•*  [w'  -  iKk'?  +  r(k';w')] 


«•.  r.  a  a 


The  new  function  r ( k ; w)  Is  related  to  G(k;w)  by  the  equation 


G(k; w)  =  - “ -  .  (30) 

[w  -  iKx  +  r(k;w)] 

Eq.  (28)  can  be  written  the  equivalent  form 

I w  +  iKk2  +  r(k;w)|2SH(k;w)  =  £(b.k'  x  k"/B)2  f(k',k") 

00 

X  |  SK(k*;w')SH(k.";w  -  w' )  dw'  .  (31) 

— 00 

which  isolates  the  nonlinear  damping  function  r(k;w).  The  WCA  equations 
are  usually  presented  in  this  form  [Eqs.  (18)  and  (19)  in  Sudan  and 
Kesklnen.  (1979)]. 

In  the  general  forms  presented  here,  neither  the  DIA  nor  the  WCA 
equations  can  be  solved,  but  we  can  establish  an  important  property  of 
Eqs.  (22)  and  (28).  Working  from  Eq.  (28),  we  define  the  argument  of 
the  summation  as  R(k,  k' .  k" ) ,  which  has  the  property 

R(k,k',k")  +  R(k',k,k")  +  R(k,k",k')  =  0  .  (32) 

It  follows  from  Eq.  (32)  that  the  summation  of  the  right-hand  side  of 
Eq.  (28)  over  all  k,  allowing  for  the  triangle  equality  constraint,  is 
zero.  Thus,  the  conservation  property  of  the  basic  equations  is 
preserved  in  the  DIA  or  WCA  equations  for  the  spectral  density  function 
itself.  This  property  is  used  to  determine  the  spectral 
characteristics . 
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SECTION  3 

DETERMINATION  OF  SDF  CHARACTERISTICS 


Sudan  and  Kesklnen  [1979]  have  argued  that  SN(k;w)  will  be  highly 
peaked  about  w  =  0  with  a  characteristic  width  that  depends  on  r(k;w). 
They  have  postulated  a  functional  relation  of  the  form 

S(k;w)  =  ♦(k)g(w/rfc)  ,  (33) 

where  Is  the  direction  averaged  value  of  r(k;w).  For  the  moment,  let 
us  consider  only  the  fact  that  we  have  effectively  assumed  that  the  time 
or  temporal  frequency  dependence  of  S  (k;w)  is  a  multiplicative  func- 

H 

tlon.  With  this  assumption  alone,  we  can  rewrite  Eq.  (28)  as 


311*1 

at 


+  Kk2*(k)  =  20(k)  23(b  X  k' 


x  k"/B)2f(k';k")4(k')4(k") 


-  4(k)  r(k) 


(34) 


where 


0(k)  =2»1 


JJ 


g[wVrfcg(w  -  w)/rfc] 
[w  +  iKk2  +  r(k;w)]* 


dw  dw' 

2t  2t 


(35) 


and  r(k)  denotes  the  value  of  -ir(k,w)  integrated  over  all  w.  Both  8(k) 
and  r(k)  are  purely  real. 

At  this  point,  the  usual  approach  is  to  assume  that  4(11)  is  Inde¬ 
pendent  of  the  direction  of  k,  as  would  be  strictly  true  when  Eq.  (26) 
holds.  Indeed,  this  same  property  would  then  apply  to  r(k)  and  6(k). 

The  remaining  term  in  Eq.  (34)  can  be  evaluated  and  a  conservation 


relation  formulated  In  terms  of  the  magnitude  k.  One  then  characterizes 
the  transfer  of  scalar  variance  from  large  to  small  wave  numbers.  It 
happens,  however,  that  the  second  term  in  Eq.  (34)  dominates  this 
process,  and  a  differential  equation  of  the  form 

k  k  k 

je(k' )dk'  +  |  r(k')4(k*)dk*  =  J  k'2K  4(k')dk'  .  (36) 

o  ak  o 

where  a  is  a  small  positive  number  results.  The  lower  limit  of  the 
nonlinear  damping  term  is  nonzero  because  arbitrarily  large  structures 
are  essentially  background  lnhomogenitles. 

In  the  stationary  state,  the  time  derivative  is  replaced  by  a 
constant.  If  K  =  0,  it  follows  that  r(k)4(k)  must  be  a  function,  say 
Y(k) ,  with  the  property 


Y(k)  *  aY(ak)  ,  (37) 

so  that  the  integral  over  ak  to  k  is  zero.  The  solution  to  Eq.  (37)  Is 
proportional  to  k  .  To  determine  the  form  of  #(k),  we  must  establish  a 
relation  between  r(fc)  and  #(k).  Sudan  and  Kesklnen  [1979]  have  derived 
an  approximate  form  by  using  the  Gaussian  form 

g(w/rj  =  [l/(2tr2)1/2|exp{-w2/r2}  (38) 

for  g(w/rfc).  With  appropriate  assumptions,  they  were  able  to  derive  the 
approximate  relation  (for  two  dimensions) 

r,  ■  v\1/2  •  (3,) 

where  I  denotes  the  angle  averaged  value  of  ♦  (k).  It  follows  that  I 
r  ^  2  R 

is  proportional  to  k  .  This  would  seem  to  imply  that  in  the  "cascade" 


range,  the  one-dimensional  spectral  density  function  is  proportional  to 
k  \  which  seems  to  be  close  to  the  measured  value.  As  we  have  noted. 


however,  this  approach  Is  Inconsistent  with  the  anisotropy  of  the 
collisional  Instability. 

We  suggest  applying  the  argument  in  a  different  way.  Because  the 
problem  is  approximately  one  dimensional,  with  most  of  the  small-scale 
structure  characterized  by  Fourier  modes  along  D,  It  seems  reasonable  to 
replace  the  angle  averaged  quantities  by  integrals  over  the  direction 
orthogonal  to  D,  which  we  take  to  lie  along  the  x  axis.  We  obtain  a 
relation  identical  to  Eq.  (36)  with  ®(k)  replaced  by  e(ky).  It  follows 
that  in  the  cascade  range, 

e(*y)r(*y]  *  *y  (40) 

and 

^y)  =  V*y/2«1/2(*y)  (41) 

from  which  it  follows  that  ®(ky)  =  fc^1;  however,  because  of  the  one- 
dimensional  approximation,  this  result  should  be  regarded  as  a  limiting 
form. 
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SECTION  4 


DISCUSSION 

We  have  shown  that  a  consistent  application  of  statistical  tur¬ 
bulence  theory  to  convective  plasma  instabilities  must  accommodate  the 
anisotropy  of  the  structure,  which  is  unique  to  the  colllsional  limit. 

To  get  some  handle  on  how  this  anisotropy  might  affect  the  theory,  we 
have  followed  the  general  prescription  suggested  by  Sudan,  but  attempted 
to  relax  the  isotropy  assumption  by  going  to  the  opposite  extreme  of  a 
one-dimensional  transfer  of  scalar  variance.  For  well-developed 
structures  transverse  to  the  dominant  convection  direction,  the  limiting 
form  of  the  one-dimensional  spectral  density  function  is  k  1.  Along  the 

convection  direction,  the  one-dimensional  spectrum  should  be  "shape- 

-2 

dominated"  and,  therefore,  k  at  long  wavelengths,  giving  way  to  a  much 
steeper  power  law  wnere  diffusion  is  effective.  Indeed,  in  the 
diffusion  range  the  spectrum  should  be  Isotropic,  at  least  if  classical 
cross-field  diffusion  is  the  dominant  mechanism  for  structure  removal. 

Ott  and  Farlev  [1974]  first  pointed  out  that  a  k  1  one-dimensional 

spectral  density  function  for  density  structure  is  the  form  that  results 

from  the  purely  dimensional  arguments  of  turbulence  theory;  however, 

they  questioned  the  existence  of  a  truly  lnertialllke  subrange  because 

it  implied  no  dependence  on  the  initial  gradient;  moreover,  the  avail- 

-2 

able  data  then  seemed  to  imply  a  k  form.  In  our  development,  we  have 
not  addressed  the  question  of  the  constants  or  other  factors  that  deter¬ 
mine  the  level  of  the  spectrum  in  addition  to  its  shape.  The  dependence 
on  a  length  scale  comes  from  the  a  in  Eq.  (36).  For  fluid  turbulence,  a 
is  related  to  the  Kolmogoroff  constant,  but  the  theory  as  it  is  applied 
here  does  not  determine  its  value  nor  does  it  identify  an  actual  physi¬ 
cal  scale  that  a  can  be  identified  with,  although  an  outer  scale  or 
maximum  size  for  transfer  of  scalar  variance  is  reasonable. 
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We  also  note  that  Batchelor  [1959],  In  treating  the  structure  of 
passive  scalar  additives  in  turbulent  fluids,  showed  that  in  the  pres¬ 
ence  of  a  linear  velocity  shear  the  average  one-dimensional  spectrum  for 
a  conserved  scalar  governed  by  Eq.  (1)  is  also  k  1.  In  fluid  turbu¬ 
lence,  this  comes  about  when  the  viscous  cutoff  for  velocity  turbulence 
occurs  at  a  larger  scale  than  the  diffusive  cutoff  for  the  scalar 
additive.  This  is  referred  to  as  the  viscous  convective  subrange 
fTennekes  and  Lumelv.  1972].  The  velocity  shear  is  a  consequence  of  the 
turbulent  dissipation  mechanism.  To  apply  Batchelor's  theory  to  convec¬ 
tive  Instabilities,  the  presence  of  the  velocity  shears  must  be  estab¬ 
lished  Independently.  The  direct  application  of  the  turbulence  theory 
is  more  satisfying. 

There  is  yet  another  turbulence  theory  argument  that  gives  rise  to 
limiting  spectral  forms  that  have  power-law  slopes  less  than  two.  These 
are  discussed  in  Kintner  and  Sevier  [1985],  and  they  involve  a  dual 
cascade,  which  can  occur  in  two  dimensional  fluids.  Indeed,  the  fact 
that  Eq.  (4)  conserves  enstropy  introduces  a  second  conserved  quantity 
because  energy  is  conserved  as  well.  Kralchnan  [1967]  showed  that  under 
these  conditions  two  types  of  cascade  can  occur.  At  scale  sizes  above  a 
stirring  scale  where  the  structure  is  Initiated,  there  is  an  upward 
cascade  of  enstropy  and  scaling  arguments  give  a  k  1  form  for  the  enerpv 
spectrum.  Below  the  stirring  scale,  there  is  a  downward  cascade  of 
energy  and  the  conventional  h  5^3  spectrum  results.  To  the  extent  that 
density  is  purely  a  passive  scalar  in  this  process  and  is  related  to 
velocity  via  Eqs.  (21)  and  (27),  its  spectrum  assumes  the  same  two- 
component  spectral  form.  The  corresponding  forms  that  apply  when  Eqs. 
(21)  and  (26),  which  are  exact,  hold  are  given  in  Table  I  of  Kintner  and 
Sevier.  [1985].  The  problem  here  Is  that  the  conditions  under  which  Eq. 
(27)  can  be  used  give  rise  to  anisotropy,  in  any  case,  would  preclude 
the  possibility  of  density  behaving  as  a  passive  scalar  as  well. 

Indeed.  Sudans's  theory  has  yet  to  be  generalized  to  accommodate  more 
than  a  single  conserved  quantity.  More  important,  the  experimental 
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data,  although  showing  two-component  forms  with  fairly  shallow  low- 
frequency  indices,  have  not  shown  the  two-component  form  predicted  by 
the  dual  cascade  theory. 

To  summarize,  the  two-dimensional  spectral  density  function  for 
convective  Instabilities  as  modeled  by  Eqs.  (1)  and  (2)  is  anisotropic 
at  the  longest  wavelengths  with  a  power-law  index  approaching  -1  in  the 
direction  perpendicular  to  the  convection  direction  and  steeping  to  -2 
or  greater  along  the  convection  direction.  Beyond  the  freezing  scale 
where  diffusion  becomes  effective,  the  spectrum  would  be  more  nearly 
Isotropic.  The  ramifications  of  this  model  for  propagation  models  are 
probably  not  severe  because  the  slope  change  occurs  only  over  a  limited 
scale-size  regime.  As  a  means  of  testing  this  hypothesis,  however,  the 
propagation  Issues  should  be  pursued. 
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ATTN.  SLCAS-AE-E 

U  S  ARMY  COMMUNICATIONS  R&D  COMMAND 
ATTN:  DRDCO-COM-RY  W  KESSELMAN 

U  S  ARMY  FOREIGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD 

U  S  ARMY  MATERIAL  COMMAND 
ATTN:  DRCLDC  J  BENDER 

U  S  ARMY  MISSILE  COMMAND 

ATTN:  DRSMI-YSO  J  GAMBLE 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN:  LIBRARY 

U  S  ARMY  SATELLITE  COMM  AGENCY 
ATTN;  AMCPM-SC  3 
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U  S  ARMY  STRATEGIC  DEFENSE  CMD 
ATTN:  DASD-H-SAV 

U  S  ARMY  STRATEGIC  DEFENSE  COMMAND 
ATTN:  ATC-OW  DAVIES 
ATTN:  ATC-R  W  DICKSON 

U  S  ARMY  TRADOC  SYS  ANALYSIS  ACTVY 
ATTN:  ATAA-PL 
ATTN:  ATAA-TCC  F  PAYAN  JR 
ATTN:  ATAA  TDC 

US  ARMY  WHITE  SANDS  MISSILE  RANGE 
ATTN:  STEWS-TE-N  K  CUMMINGS 

DEPARTMENT  OF  THE  NAVY 

COMMAND  &  CONTROL  PROGRAMS 
ATTN:  OP  941 

JOINT  CRUISE  MISSILES  PROJECT  OFC  (PM  3) 
ATTN:  JCMG  707 

NAVAL  AIR  SYSTEMS  COMMAND 
ATTN.  PMA271 

NAVAL  ELECTRONICS  ENGRG  ACTVY,  PACIFIC 
ATTN:  CODE  250  D  OBRYHIM 

NAVAL  INTELLIGENCE  SUPPORT  CTR 
ATTN:  NISC-50 

NAVAL  OCEAN  SYSTEMS  CENTER 
ATTN:  CODE  532 
ATTN:  CODE  54  J  FERGUSON 

NAVAL  RESEARCH  LABORATORY 

ATTN:  CODE  4180  J  GOODMAN 
ATTN:  CODE  4700  S  OSSAKOW 
ATTN:  CODI.  4720  J  DAVIS 
ATTN:  CODE  4732  BRIPIN 
ATTN:  CODE  4750  P  RODRIGUEZ 
ATTN.  CODE  4780  J  HUBA 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN  CODE  H  21 

NAVAL  UNDERWATER  SYS  CENTER 
ATTN:  CODE  3411,  J  KATAN 

OFC  OF  THE  DEPUTY  CHIEF  OF  NAVAL  OPS 
ATTN:  NOP 94 ID 

ATTN:  OP  654(STRAT  EVAL  4  ANAL  BR) 
ATTN:  OP981N 

SPACE  4  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN:  CODE  3101  T  HUGHES 
ATTN:  CODE  501A 
ATTN:  PD50TD 
ATTN:  PDE  110  X1  B  KRUGER 


ATTN:  PDE-110-11021  G  BRUNHART 
ATTN:  PME  106-4  S  KEARNEY 
ATTN:  PME-106  F  W  DIEDERICH 

STRATEGIC  SYSTEMS  PROGRAM  OFFICE  (PM  1) 
ATTN:  NSP-L63  (TECH  LIB) 

ATTN:  NSP-2141 
ATTN:  NSP-2722 

THEATER  NUCLEAR  WARFARE  PROGRAM  OFC 
ATTN:  PMS-42331F(D  SMITH) 

DEPARTMENT  OF  THE  AIR  FORCE 

AIR  FORCE  CTR  FOR  STUDIES  4  ANALYSIS 
ATTN:  AFCSA/SAMI  (R  GRIFFIN) 

ATTN:  AFCSA/SASC 

AIR  FORCE  GEOPHYSICS  LABORATORY 
ATTN:  LID/J  RAMUSSEN 
ATTN:  LIS  J  BUCHAU 
ATTN:  LS 

ATTN:  LS/R  O’NIEL 
ATTN:  LSI/ H  GARDINER 
ATTN:  LYD/K  CHAMPION 

AIR  FORCE  SATELLITE  CTRL  FACILITY 
ATTN:  WE 

AIR  FORCE  SPACE  DIVISION 
ATTN.  YA 
ATTN:  YG 
ATTN:  YK 
2  CYS  ATTN:  YN 

AIR  FORCE  SPACE  TECHNOLOGY  CENTER 
ATTN:  XP 

AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  TN 

AIR  FORCE  WEAPONS  LABORATORY,  NTAAB 
ATTN:  NTN 
ATTN:  SUL 

AIR  FORCE  WRIGHT  AERONAUTICAL  LAB/AAAD 
ATTN:  WHUNT 

AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL-LSE 

BALLISTIC  MISSILE  OFFICE/DAA 
ATTN:  ENSN 

DEPUTY  CHIEF  OF  STAFF/AF-RDQM 
ATTN:  AF/RDQI 

DEPUTY  CHIEF  OF  STAFF/AFRDS 

ATTN  AFRDS  (SPACE  SYS  4  C3  DIR) 
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ELECTRONIC  SYSTEMS  DIVISION/SC 
ATTN:  SCS-1E 

ATTN.  SCS-2  LTC  ALEXANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
ATTN:  OCDS  R  SCHNEIBLE 
ATTN:  TSLD 

SPACE  COMMAND/DC 
ATTN:  XP 

STRATEGIC  AIR  COMMAND/NRI-STINFO 
ATTN:  NRI/STINFO 

STRATEGIC  AIR  COMMAND/XPFC 
ATTN:  XPFS 

STRATEGIC  AIR  COMMAND/XPQ 
ATTN:  XPQ 

DEPARTMENT  OF  ENERGY 

DEPARTMENT  OF  ENERGY 
ATTN:  DP-233 

EG&G,  INC 

ATTN:  D  WRIGHT 

UNIVERSITY  OF  CALIFORNIA 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  L-31  R  HAGER 
ATTN:  L  53  TECH  INFO  DEPT.  LIBRARY 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN.  D  SAPPENFIELD 
ATTN:  D  SIMONS 
ATTN:  J  WOLCOTT 
ATTN:  MSJZINN 
ATTN:  R  JEFFRIES 
ATTN:  RW  WHITAKER 
ATTN:  TKUNKLE 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  D  HARTLEY  8300 

SANDIA  NATIONAL  LABORATORIES 

ATTN.  A  D  THORNBROUGH  400 

ATTN:  D  DAHLGREN  6440 

ATTN:  ORG  1231  T  P  WRIGHT 

ATTN:  ORG  314  W  0  BROWN 

ATTN:  ORG  332  R  C  BACKSTROM 

ATTN:  SPACE  PROJECT  DIV 

ATTN:  TECH  LIB  3141  (RPTS  RECEIVING  CLRK) 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  OSWR/NED 

ATTN:  OSWR/SSD  FOR  K  FEUERPFETL 
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U  S  DEPARTMENT  OF  COMMERCE 

ATTN:  A  JEAN  (UNCL  ONLY) 

ATTN:  WUTLAUT 

US  DEPARTMENT  OF  STATE 
ATTN:  PM/TMP 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  D  OLSEN 
ATTN.  E  RODRIGUEZ 
ATTN:  IGARFUNKEL 
ATTN:  JKLUCK 
ATTN:  J STRAUS 
ATTN:  KSCHO 
ATTN.  R  SLAUGHTER 
ATTN:  T  SALMI 

AEROSPACE  CORP 

ATTN:  SMEWATERS 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
ATTN:  SECURITY 

AUSTIN  RESEARCH  ASSOCIATES 
ATTN:  J  THOMPSON 

AUTOMETRIC,  INC 

ATTN:  C LUCAS 

BDM  CORP 

ATTN:  AVITELLO 
ATTN:  L  JACOBS 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  CPRETTIE 
ATTN:  J  WORKMAN 
ATTN:  S  BRECHT 

BOEING  CO 

ATTN:  GHALL 
ATTN:  D  CLAUSON 

BR  COMMUNICATIONS 

ATTN:  J  MCLAUGHLIN 

CALIFORNIA  RESEARCH  &  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 

CHARLES  STARK  DRAPER  LAB,  INC 
ATTN:  ATETEWSKI 

COMMUNICATIONS  SATELLITE  CORP 
ATTN:  G  HYDE 

COMPUTER  SCIENCES  CORP 
ATTN:  F  EISENBARTH 
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CORNELL  UNIVERSITY 

ATTN:  D  FARLEY  JR 
ATTN:  M  KELLY 

ELECTROSPACE  SYSTEMS,  INC 
ATTN:  HLOGSTON 
ATTN:  P  PHILLIPS 

EOS  TECHNOLOGIES,  INC 
ATTN:  B  GABBARD 
ATTN:  WLELEVIER 

GENERAL  ELECTRIC  CO 

ATTN:  ASTEINMAYER 
ATTN:  CZIERDT 

GEO  CENTERS,  INC 

ATTN:  E  MARRAM 

GTE  GOVERNMENT  SYSTEMS  CORPORATION 
ATTN:  W  I  THOMPSON,  III 

HARRIS  CORP 

ATTN:  EKNICK 

HSS,  INC 

ATTN:  D HANSEN 

IBM  CORP 

ATTN:  HULANDER 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  E  BAUER 
ATTN:  HWOLFHARD 

JAYCOR 

ATTN:  J  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  CMENG 
ATTN:  JD  PHILLIPS 
ATTN:  J  NEWLAND 
ATTN  K  POTOCKI 
ATTN:  R  STOKES 
ATTN.  T  EVANS 

KAMAN  SCIENCES  CORP 
ATTN.  E  CONRAD 

KAMAN  TEMPO 

ATTN:  BGAMBILL 
ATTN:  DASIAC 
ATTN.  R RUTHERFORD 
ATTN:  W MCNAMARA 

KAMAN  TEMPO 

ATTN:  DASIAC 


LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  JKUMER 
ATTN:  R  SEARS 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
2  CYS  ATTN:  D  CHURCHILL 

M  I  T  LINCOLN  LAB 

ATTN:  D  TOWLE 
ATTN:  IKUPIEC 

MAXIM  TECHNOLOGIES.  INC 
ATTN:  J  LEHMAN 
ATTN:  J MARSHALL 
ATTN:  JSO 
ATTN:  NCIANOS 

MCDONNELL  DOUGLAS  CORP 
ATTN:  RHALPRIN 

METEOR  COMMUNICATIONS  CORP 
ATTN:  R  LEADER 

MISSION  RESEARCH  CORP 
ATTN:  BR  MILNER 
ATTN:  CLAUER 
ATTN:  CRINO 
ATTN:  D  ARCHER 
ATTN:  DKNEPP 
ATTN:  FFAJEN 
ATTN:  F  GUIGLIANO 
ATTN:  G  MCCARTOR 
ATTN:  KCOSNER 
ATTN:  RBIGONI 
ATTN.  R  BOGUSCH 
ATTN:  R  DANA 
ATTN:  R  HENDRICK 
ATTN:  RKILB 
ATTN:  R  MORGANSTERN 
ATTN:  S  GUTSCHE 
ATTN:  TECH  LIBRARY 

MITRE  CORPORATION 

ATTN:  AKYMMEL 
ATTN:  C  CALLAHAN 
ATTN:  D  RAMPTON,  PH  D 
ATTN:  DR  D  RAMPTON 
ATTN.  MRDRESP 
ATTN:  RDRESP 

MITRE  CORPORATION 

ATTN.  J  WHEELER 
ATTN:  MHORROCKS 
ATTN:  RCPESCI 
ATTN:  W  FOSTER 

NORTHWEST  RESEARCH  ASSOC,  INC 
ATTN  EFREMOUW 
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PACIFIC  SIERRA  RESEARCH  CORP 
ATTN:  E  FIELD  JR 
ATTN:  F  THOMAS 
ATTN:  H  BRODE,  CHAIRMAN  SAGE 

PHOTOMETRICS,  INC 

ATTN:  ILKOFSKY 

PHYSICAL  RESEARCH  INC 
ATTN:  H  FRITZ 
ATTN:  JJORDANO 

PHYSICAL  RESEARCH,  INC 
ATTN:  RDELIBERIS 
ATTN:  T  STEPHENS 

PHYSICAL  RESEARCH,  INC 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  W  SCHLUETER 

R  &  D  ASSOCIATES 

ATTN:  BLAMB 
ATTN:  C  GREIFINGER 
ATTN:  F  GILMORE 
ATTN:  GHOYT 
ATTN:  HORY 
ATTN:  MGANTSWEG 
ATTN:  M  GROVER 
ATTN:  RTURCO 
ATTN  WKARZAS 

R  4  D  ASSOCIATES 

ATTN:  B  WEBSTER 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DTELAGE 
ATTN:  M  CROSS 

SRI  INTERNATIONAL 

ATTN:  D  MCDANIEL  (UNCL  ONLY) 

ATTN:  WCHESNUT 
ATTN:  WJAYE 

STEWART  RADIANCE  LABORATORY 
ATTN:  RHUPPI 

TECHNOLOGY  INTERNATIONAL  CORP 
ATTN,  W  BOQUIST 

TELECOMMUNICATION  SCIENCE  ASSOCIATES 
ATTN:  R  BUCKNER 

TOYON  RESEARCH  CORP 

ATTN:  J  GARBARINO 
ATTN:  JISE 

TRW  ELECTRONICS  4  DEFENSE  SECTOR 
ATTN:  DRDGRYBOS 
ATTN:  RPLEBUCH 

UTAH  STATE  UNIVERSITY 
ATTN.  A  STEED 
ATTN:  D  BURT 

ATTN:  K  BAKER.  DIR  ATMOS  4  SPACE  SCI 
ATTN:  L  JENSEN,  ELEC  ENG  DEPT 

VISIDYNE.  INC 

ATTN  J CARPENTER 
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R  4  D  ASSOCIATES 

ATTN:  GGANONG 

RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  E  BEDROZIAN 
ATTN:  P  DAVIS 


RAND  CORP 

ATTN:  B  BENNETT 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  S  QUILICI 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  C  SMITH 
ATTN:  D  HAMLIN 
ATTN:  ESTRAKER 
ATTN:  LLINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN  R  LEAOABRAND 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN  J COCKAYNE 
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